Crossover properties from random percolation to frustrated percolation 
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We investigate the crossover properties of the frustrated percolation model on a two-dimensional 
square lattice, with asymmetric distribution of ferromagnetic and antiferromagnetic interactions. 
We determine the critical exponents v, 7, and j3 of the percolation transition of the model, for 
various values of the density of antiferromagnetic interactions tt in the range < tt < 0.5. Our 
data are consistent with the existence of a crossover from random percolation behavior for tt = 0, 
to frustrated percolation behavior, characterized by the critical exponents of the ferromagnetic 1/2- 
state Potts model, as soon as tt > 0. 
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I. INTRODUCTION 

The cluster approach was introduced by Kasteleyn and 
Fortuin (KF) fjj and Coniglio and Klein (CK) § in ferro- 
magnetic spin systems. In the Coniglio-Klein approach, 
one puts a bond between two nearest neighbor (NN) spins 
with probability p = 1 — e~ 2,3J if the spins are parallel, 
where J is the spin interaction and f3 = l/fc^T, and 
with probability zero if they are antiparallel. In this way 
a bond configuration C on the lattice has a statistical 
weight 



W(C) = e ^q N(C) , 



(1) 



where \x — log (i^jjl — log(e 2 ^ J — 1) is the chemical po- 
tential of the bonds, b(C) is the number of bonds and 
N(C) the number of clusters of the configuration C, and 
q is the multiplicity of the spins (q = 2 for Ising spins). 
This defines a percolation model, in which clusters per- 
colate at the ferromagnetic critical point, with the same 
critical indices of the original spin model. 

The KF-CK approach has been extended also to frus- 
trated spin models, as for example spin glasses ||. In 
these models the disorder is produced by a quenched dis- 
tribution of antiferromagnetic (—J) and ferromagnetic 
(J) interactions on the lattice. Just like in the ferromag- 
netic case, one puts a bond between two NN spins if they 
satisfy the interaction, with a probability p = 1 — e~ 2/3J . 
The main difference here is that, due to frustration, the 
spins cannot satisfy simultaneously all the interactions 
on the lattice. More specifically, if a closed path on the 
lattice contains an odd number of antiferromagnetic in- 
teractions, not all the interactions belonging to the path 
can be satisfied simultaneously. Such a path is called 
"frustrated loop" , and since bonds can be put only be- 
tween spins that satisfy the interactions, a frustrated loop 
cannot be completely occupied by bonds. Therefore, the 
statistical weight of a bond configuration C will be now 



W(C) 



e fj.b(C) q N(C) 



if C is not frustrated, 
if C is frustrated, 



(2) 



where C is said frustrated if it contains one or more frus- 
trated loops completely occupied by bonds. 

It has been shown by renormalization group meth- 
ods on a hierarchical lattice Q that this model exhibits 
two phase transitions, for every value of the multiplic- 
ity q of the spins. The first transition, at a tempera- 
ture Tsa{q)i is in the universality class of the Ising SG 
transition, while the other transition, at a temperature 
T P (q) > Tsciq), is a percolation transition in the univer- 
sality class of the ferromagnetic q/2-state Potts model 
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The frustrated percolation model has proven to be a 
suitable model to the study of complex systems, such as 
spin glasses Q , glasses Q and in general all those systems 
in which connectivity and frustration play a fundamental 
role, (for a review see ||). 

The aim of the present paper is the study of the perco- 
lation transition of the model with q — 1, the "bond frus- 
trated percolation model" , for a variable density tt of an- 
tiferromagnetic interactions in the interval < tt < 0.5. 
We determine the critical probability p c {tt) and the criti- 
cal exponents v(tt), (3(tt), and 7(71"), by performing Monte 
Carlo simulations on lattices of different size, and us- 
ing scaling laws (see Sect. III). Finally the results are 
compared with the theoretical predictions for the two 
extreme cases, the pure ferromagnetic case tt = 0, that 
corresponds to random bond percolation (1/f = 0.75, 
f3/v = 0.1042, j/i/ = 1.7917), and the symmetric case 
tt = 0.5, that corresponds to the 1/2-state Potts model 
(1/u = 0.5611, (3/v = 0.08276, 7/1/ = 1.8346). 



II. DEFINITION OF THE FRUSTRATED 
PERCOLATION MODEL 

The bond frustrated percolation model can be defined 
in the following way. Consider a two-dimensional square 
lattice, with NN interactions between sites. These inter- 
actions can be ferromagnetic with probability 1 — tt and 
antiferromagnetic with probability tt. The distribution of 
interactions is quenched, so it is set at the beginning and 
does not evolve with time in the dynamics of the system. 
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Each edge of the lattice, connecting a pair of NN sites, 
can be connected by a bond or not, and the state of the 
system is completely specified by the bond configuration. 
We give to a bond configuration C a statistical weight 



W(C) = 



eM fc(C) if C is not f rus trated, 
if C is frustrated, 



(3) 



where fi — log ( jz^j , p is a probability that is connected 

to the temperature via the relation p = 1 — e~ 2,3J , and 
6(C) is the number of bonds of the configuration C. 

The presence of frustration induces a complex behav- 
ior in both the static and dynamic properties of the sys- 
tem, with the presence of many metastable states, and 
high free energy barriers separating them. For tt — the 
model coincides with random percolation. For tt = 1/2, 
that is equal density of ferromagnetic and antiferromag- 
netic interactions, the system is expected to have two 
phase transitions Q. The first at lower temperature in 
the same universality class of the Ising spin glass tran- 
sition, which in two dimensions is at T = 0, that is at 
p = 1. Below this temperature the system is frozen, er- 
godicity is broken and the system remains trapped in a 
finite region of phase space. The second transition at a 
higher temperature (lower probability) , is the percolation 
transition of the cluster of bonds, and belongs to a differ- 
ent universality class, namely that of the ferromagnetic 
1/2-state Potts model. 

We will study the percolation transition as a function 
of the density of antiferromagnetic interactions tt. We 
will see that a very small amount of antiferromagnetic 
interactions is already sufficient to change the universal- 
ity class of the transition. 

III. MONTE CARLO RESULTS 

We have studied the percolation transition of the 
model on a 2D square lattice with periodic boundary 
conditions, for different values of tt (0 < tt < 0.5) and 
for different lattice sizes L (L — 16,24,32,40, occasion- 
ally larger). For each size of the lattice a configuration 
of interactions is preliminary produced by setting ran- 
domly on the lattice a number N a of antiferromagnetic 
interactions and Nf = 2L 2 — N a of ferromagnetic inter- 
actions. The corresponding density of antiferromagnetic 
interactions is given by tt = N a /2L 2 . 

We have determined the critical probability of perco- 
lation Pc(tt) and the critical exponents f(7r), 7(71"), and 
/3(7r), by performing Monte Carlo simulations and using 
scaling laws. In order to generate bonds configurations 
with the appropriate weights we used an algorithm de- 
vised by Sweeny p|, s uitably modified to treat the frus- 
tration occurrence ]10| . For each tt the quantities of inter- 
est were averaged on many configurations of interactions. 

For each value of tt the critical probability of perco- 
lation has been determined using the following scaling 



law for the probability of a spanning cluster being 
present on the lattice |ll[ 

P oo (L,p)=P 0O [L l ^(p-p c )], lorL^oo,p^p c . (4) 

Therefore p c is found to be the point where the curves of 
Pqo &S cL function of p, with different lattice sizes, inter- 
sect. 

In Fig. n are shown the curves P^L^p) for tt = 0.5, 
and in Fig. 2j the values obtained for p c (tt). The intersec- 
tion point is clearly singled out. This procedure enabled 
us to determine p c with an indetermination of ±0.001, or 
in the worst cases of ±0.002. 

The simulations were performed on lattices of size 
L = 16, 24, 32, 40, for each L about 500 MC steps 
were produced to thermalize the system whereas 20, 000- 
30, 000 MC steps were used to average. Moreover for each 
value of tt we have averaged on a number n of configura- 
tions of interactions ranging from n = 80 for L = 16 to 
n = 30 for L = 40. 

The scaling law in Eq. (|4|) enables us to get the value 
of the exponent 1 jv as well, once p c has been determined, 
by choosing the value which gives the best data collapse 
of the curves (see Fig. |3|). In Fig. |^ we give the values 
of 1/v(tt) obtained. As tt increases from zero to positive 
values, a sudden change in the universality class of the 
transition can be seen. A crossover region extending from 
tt = to tt ~ 0.05 occurs. From this point the value of v 
can be considered in agreement with the predicted value 
for the 1/2-state Potts model {1/v = 0.56, marked by 
the horizontal straight line). 

The errors A(l/i>) and Ap c were computed as the am- 
plitudes of the intervals in the values of \jv and p c for 
which a good data collapse was obtained. We remark 
that these errors do not take into account the finiteness 
of the system. Thus the values for the critical exponents 
must be regarded as effective values which would give 
correct results only in the asyntotic limit (L — » 00). 

The mean cluster size x is defined as 



V ^ 



s 2 n„ 



(5) 



where V = L 2 is the number of sites on the lattice, n s 
is the number of clusters of size s, and the sum extends 
over the cluster sizes s. In the thermodynamic limit L — > 
00 the mean cluster size diverges as \p — p c | -7 , when 
the probability p approaches its critical value. For finite 
systems, \ obeys a finite size scaling Jll|] 

X (L,p) = L-y^x[L 1/u (P'Pc)} forL^co,p^ Pc . (6) 

The density of the largest cluster plays the role of the 
order parameter in the system, being zero for p < p c in 
the thermodynamic limit, and oc {p — p c )^ for p > p c - 
This quantity as well obeys a finite size scaling 



Poo(L,p) = L-^p^L^ip-pc)} for L > 00. 



P->Pc 



(7) 
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Therefore simulating the system at the computed value 
of p c enables us to get 7 and /3 from a log-log plot of the 
relations x{L,p = p c ) oc L 1 ^ , Poo(L,p = p c ) oc L~^l v . 
In Fig. ^|one such plot is shown, for tt = 0, while in Fig. ^ 
and Fig. [?] are shown the plots of 7/1/ and (3/v obtained, 
for different values of tt. The horizontal straight lines 
mark the predicted values of the exponents, 7/z/ = 1.83 
and (3/v = 0.083, corresponding to the ferromagnetic 
1/2-state Potts model. We see that for tt > 0.05 the 
computed value of the exponents is in good agreement 
with the prediction. 

These results are consistent with the picture that there 
are two universality classes: random percolation at tt = 
0, and frustrated percolation for tt > 0. Data at our 
disposal cannot exclude however the possibility that the 
random percolation behavior extends from tt = to a 
value tt* smaller than 0.05. 



IV. CONCLUSIONS 

We have investigated the percolation transition of the 
asymmetric frustrated percolation model in two dimen- 
sions by using Monte Carlo simulation. 

From the analysis of critical exponents 7(71"), 
(3(tt), in the interval < tt < 0.5, it seems reasonable 
to assume that a very small concentration (tt ~ 0.05) of 
antifcrromagnetic interactions is already sufficient to pro- 
duce the change in the universality class of the transition. 
This means that the effects of disorder and frustration are 
important even for such small values of tt. 

Our results are consistent with the existence of a sharp 
crossover from random percolation for tt = 0, to frus- 
trated percolation, characterized by the exponents of the 
ferromagnetic 1/2-state Potts model, as soon as tt > 0. 
However we cannot rule out numerically the presence of 
a tricritical point, at low values of tt, dividing random 
percolation from frustrated percolation exponents. 
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FIG. 1. Probability of percolation Poo, as a function of the 
probability p, for lattice sizes L = 16, 24, 40 and for n — 0.5. 
The critical probability p c is found to be the intersection point 
of the different curves. 
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FIG. 2. The critical probability of percolation p c as a func- 
tion of the density of antiferromagnetic interactions tt. The 
value p c = 0.5 found for n = is in agreement with the value 
of the random percolation. 
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FIG. 3. Scaling behavior of the probability of percolation 
Poo, for 7T = 0.5. The exponent 1/v is found as the value 
which gives the best data collapse. The value 1/u = 0.56 
found is in agreement with the value of the 1/2-state Potts 
model. 



FIG. 5. Log-log plot of the mean cluster size \ m function 
of L, at p c and for tt — 0.5, and linear fit of the data. The 
critical exponent y/v is found to be the slope of the straight 
line. 
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FIG. 4. Plot of the exponent 1/v in function of the density 
of antiferromagnetic interactions 7r. The horizontal straight 
line represents the value of 1/v for the frustrated percolation. 



FIG. 6. Plot of the exponent 7/1/ in function of the density 
of antiferromagnetic interactions 7r. The horizontal straight 
line represents the value of 7/1/ for the frustrated percolation. 
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FIG. 7. Plot of the exponent fi/v in function of the density 
of antiferromagnetic interactions ■n. The horizontal straight 
line represents the value of (3/v for the frustrated percolation. 
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